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' ABSTRACT 

We present the results of an observing program with the SCUBA bolometer array to 
measure the submillimetre (submm) dust continuum emission of 24 distant {z > 1) 
radio galaxies. We detected submm emission in 12 galaxies with S/N > 3, including 9 
detections at z > 3. When added to previous published results these data almost triple 
the number of radio galaxies with z > 3 detected in the submm and yield a sample 
of 69 observed radio galaxies over the redshift range z = 1-5. We find that the range 
. in rest-frame far-infrared luminosities is about a factor of 10. We have investigated 

' the origin of this dispersion, correlating the luminosities with radio source power, 

\ size, spectral index, K-bemd magnitude and Lya luminosity. No strong correlations 

are apparent in the combined data set. We confirm and strengthen the result from 
previous submm observations of radio galaxies that the detection rate is a strong 
function of redshift. We compare the redshift dependence of the submm properties 
of radio galaxies with those of quasars and find that for both classes of objects the 
observed submm flux density increases with redshift to z w 4, beyond which, for 
the galaxies, we find tentative evidence for a decline. We find evidence for an anti- 
correlation between submm luminosity and UV polarisation fraction, for a subsamplc 
of 13 radio galaxies, indicating that starbursts are the dominant source of heating for 
dust in radio galaxies. 
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1 INTRODUCTION 

There is strong evidence that powerful high redshift radio 
galaxies (HzRGs; z > 2) are the progenitors of the bright- 
est cluster ellipticals seen today. HzRGs are the infrared 
brighte st and presumably the most massive galaxies at any 
epoch foe Breuck et afll2002l) and host actively-accreting 
super massi ve black holes with masses of order 10 9 M© 
jLacv et alJl200 it iDunlop et al]l2003h . Therefore, they are 
key objects for studying the formation and evolution of mas- 
sive galaxies and super-massive black holes. 

HzRGs are likely to be in an important phase of 
their formation process for several reasons: They have large 
reservoirs of gas from which they could be forming, as 
show n by spectacular (> 100 kpc) luminous Lya haloes 
(e.g., lMcCarthvlll993 Ivan Oiik et aHll996t | Rcul and et alJ 
2003) and widespread H I ab sorption features in the Lya 
profiles (Ivan Oiik et al.lll997i) . Their rest-frame UV mor- 



phologies are characterized by clumpy structures, simi- 
lar to the Lyman-break galaxies at z ~ 3 , that will 
merge with the central galaxy on d ynamical time-scales of 
10 8 yrs JPentericci et al.ll998tll999lh In the case of 4C 41.17 
there is direct evidence for massive star formation (up 
to ~ 1500 Mq yr _1 after c orrection for ex tinction) based 
on stellar absorption-lines foev et al]ll997l) . Finally, mm- 
interferometry studies of CO line and continuum emission 
for three z > 3 HzRGs have shown that the star for- 
mation occurs galaxy wide over distances up to 30 kpc 
foapadopoulos et all200ol ; lDe Breuck et al .120031) and there 
is even evidence fo r star formation on scales of 250 kpc 
JStevens et aljr2 003 ) . Together this suggests that we are ob- 
serving not merely scaled up versions of local ultraluminous 
infrared galaxies (ULIRGs) where the bursts are confined to 
the inner few kpc, but wide-spread starbursts within which 
the galaxies are forming the bulk of their eventual stellar 
populations. 
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HzRGs are an important sample for studying the star 
formation history of the universe because their selection is 
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based on long wavelength radio emission whose propagation 
is not affected by the presence of dust. Dust is expected 
to play a significant role in star forming regions, absorb- 
ing UV/optical radiation from the star burst and reradiat- 
ing it at far-infrared (FIR) wavelengths iSanders fc Mirabell 
Il99firl . Optical searches for di stant galaxies (e.g., us- 
ing the Lyman-br eak technique; ISteidel et alJ Il996l . Il999t 
lOuchi et al.l 1200 J) are thus likely to be biased against 
dusty objects. Finding distant star forming galaxies through 
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most obscured sources. So far there has been little over- 
lap between the optical and submm selected star forming 
sources (selection o n very red near-IR colours may prove 
more fruitful; e.g., iFraver et aDl2004) . It remains unclear 
whether they are members of a continuous popu lation (e.g., 
lAdelberger fc Steidell Eoobl IWebb et al.l l2003bl) and argu- 
ments have been made that either co uld dominate the 
star formation density at high redshift felain et alJll999t 
lAdelberger fc SteidejExibT) . Since selection at radio wave- 
lengths circumvents the aforementioned selection biases it 
could help determine the relative contributions of obscured 
and unobscured star formation to the star formation history 
of t he universe. 

lArchibald et alJ i200ll . hereafter A01) have conducted 
the first systematic submm survey to study the star for- 
mation history of radio galaxies over a redshift interval of 
0.7 < z < 4.4. In their sample of 47 galaxies, they found 
evidence for a considerable range in FIR luminosities, a sub- 
stantial increase in 850 ^.m detection rate with redshift and 
that the average 850 urn luminosity rises at a rate (1 + z) 3-4 
out to z ~ 4. These results prompt the following questions: 
Is the dispersion in FIR luminosities due to differences in 
their star formation rates or dust contents? Does the strong 
increase with redshift reflect an increase in star formation 
rates or could it be related to changing dust properties? Does 
the FIR luminosity keep on rising with redshift or does it 
level off and is there a redshift cut-off? Are the inferred star 
formation rates comparable to those derived from the opti- 
cal/UV? Do the submm properties of quasars (QSOs) and 
radio galaxies show similar trends or do the two classes of 
objects evolve differently? 

The submm findings from A01 were based on a limited 
number of detections at high redshift (z > 3). To put these 
results on a statistically firmer footing and search for possi- 
ble correlations with other galaxy parameters more submm 
observations were required. Here we present such observa- 
tions of all z > 3 HzRGs k nown at the beginning of 2001 
(e.g., iDe Breuck et alJl200lTl which had not been observed 
in the submm. Adding these to the survey of A01 almost 
triples the number of detections at high redshift, creating a 
sample which is statistically significant over the full redshift 
range z = 1 - 5. 

The structure of this paper is a follows: the sample se- 
lection, observations and data analysis are described in Sec- 
tion 2. Results and notes on some individual sources are 
presented in Section 3. Various correlations with submm 
properties of HzRGs are investigated in Section 4 and de- 
scribed in detail in Section 5. Section 6 presents a com- 
parison between HzRGs and QSOs. We discuss and sum- 



marize our conclusions in Section 7. Throughout this pa- 
per, we adopt a flat universe with S2m = 0.3, Qa = 0.7, 
and Ho = 65 km s" 1 Mpc" 1 . Using this cosmology the look- 
back time at z ~ 2.5 (the median redshift of our sample) 
is 11.7 Gyr and a galaxy at such a redshift must be less 
than 2.8 Gyr old. 



2 SAMPLE SELECTION AND OBSERVATIONS 

The observations presented here include submm observa- 
tions of distant radio galaxies. The targets were selected 
from an increasing sample of HzRGs that i s the result of an 
ongoing effo rt by our group and others lY)e Breuck et all 
2000a, 1200 ll . de Vries et al. in preparation; Spinrad pri- 
vate communication) to find distant radio galaxies based 
on Ultra Steep Radio Spectrum (USS; i.e. red radio color) 
and near-IR id entification selection criteria (for details see 
IDe Breuck et alJl200ll) . 

We selected all HzRGs known at the beginning of 2001 
with redshifts z > 3 and declination S > —30° that did 
not have prior submm observations. Our aim was to observe 
a significant sample of HzRGs to complement the observa- 
tions of A01 and, in particular, to obtain better statistics 
at the highest redshifts. MG 2141+192 and 4C 60.07 were 
observed in both programs, because, at the time of observa- 
tion, their inclusion in the A01 sample was unknown to us. 
The 850 ^m results for B3 J2330+3927, 6C J1908+722, and 
4C 60.07 have b een published previously as part of their CO 
imag ing studies iPapadopoulos et alJl200ol : lDe Breuck et alJ 
1200311 . MRC 1138-262 was included in t he program be- 
cause of its wealth of sup porting data (e.g., iPentericci et alj 
119971: ICarilli et al]l2002T) and WN J1115+5016 because it 
is one of only two radio galaxies showing a bro ad absorb- 
tion line (BAL) system JPe Breuck et alJlioOlh . the other 
BAL radio galaxy, 6C J190 8+722, being a strong CO emitter 
JPapadopoulos et aljlioool) . WN J0528+6549 at z = 1.210 
was observed because it was first thought to be at redshift 
z — 3.120 (actually belonging to another galaxy on the slit). 

The coordinates, redshifts, largest angular sizes of the 
radio sources, Lya fluxes, /('-band magnitudes, and refer- 
ences for the full sample observed in the submm are listed 
in Table 



2.1 SCUBA photometry 

The observations were carried out between October 1997 
and January 2002 with the Submillimet re Co mmon-User 
Bolometer Array (SCUBA: [Holland et allll999l) at the 15 m 
James Clerk Maxwell Telescope (JCMT). We observed at 
450 \im and 850 [im wavelengths resulting in beam sizes of 
7.5" and 14.7" respectively. We employed the 9-point jiggle 
photometry mode, which samples a 3 x 3 grid with 2" spac- 
ing between grid points, while chopping 45" in azimuth at 
7.8 Hz. Frequent pointing checks were performed to ensure 
pointings better than 2" and reach optimal sensitivity. 

Our original goal was to observe all sources down to 
1 mjy rms at 850 urn. This is a sensible limit, since at 850 u.m 
confusion becomes a problem fo r sources weaker than 2 mjy 
jHuehes et alJll998l:lHogj200ll) . it is obtainable in 3 hrs per 
source, and it is matched to the survey by A01. However, 
because of scheduling constraints, and because our priority 
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Table 1. Rcdshifts, radio positions, largest angular sizes, Ly« fluxes, A-band magnitudes and references to papers from 
which these data were taken for all objects that were observed in our submm program. The Lya fluxes are in units of 
10 — 16 ergs -1 cm -2 , and the if -band magnitudes were measured in a 64 kpc diameter aperture where possible. References: 
DB99.DB00a.DB00b.DB01.DB02.DB03a.DB03b= lDe Breuck et all ll99Sl l2000albl I200lll200ll200l De Breuc k et al. in 
preparation), dV03 = d e Vries et al . in preparation, Dey9 9 = IPevI il999 | ). Dri97 = iDrinkwater et all Jl997f). ER96 = 
lEales fc Rawlingsl il996h. McC90 = IMcCarthv et all Il990fl. McC96 = IMcCarthv et all il996ft. Kap9 8 = iKapahi et all 
Il998f). PapOO = iPapadoDoulos et all fcOOdV P en97 = iPentericci et all Il997tl. R en98 = iRengelinkl il99Sft . Rot97 = 
R.ott gering et all Il997tl . S99 = IStern et all ll99<t) . vB99 = Ivan Breugel et all il99Sh . 
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was to obtain a large sample of HzRGs with 850 itm detec- 
tions, this limit was not always reached. Rather, the next 
target was observed as soon as an apparent 5a detection 
had been obtained at 850 (im. 

The atmospheric optical depths rsso, T450 were calcu- 
lated using the empir ical CSO-tau correlations given by 
I Archibald et al.| ll2002t) . unless the values obtained through 
skydips disagreed strongly, in which case those were used in- 
stead. The optical depth rsso varied between 0.14 and 0.38 
with an average value of 0.26. The data were clipped at the 
4t7 level to ensure accurate determination of the sky level, 
flat-fielded, corrected for extinction, sky noise was removed 
after which they were co-added and clipped at the 2.5cr 
level using t he Scuba User Reduction F acility software pack- 
age fSURF: Ijenness fc Lightfootill9 98P). following standard 
procedures outlined in the SCUBA Photometry Cookbook 1 . 
The concatenated data were checked for internal consistency 
using a Kolmogorov-Smirnov (K-S) test and severely devi- 
ating measurements (if any) were removed. Finally, flux cali- 
bration was performed using HLTAU, OH231.8 and CRL618 
as photometric calibrators. The typical photometric uncer- 
tainty for our program is of order 10-15 per cent, as esti- 
mated from our results on three sources (TN J0121+1320, 



1 The SCUBA Photometry Cookbook is available at 
http: / /www. star link. rl.ac.uk/star/docs/scl0.htx/scl0. html 



TN J1338-1942, and MG 2141+192) that were observed 
at two separate instances each. This photometric accu- 
racy is consistent with an estimated 10 per cent system- 
atic uncertainty in the 850 u.m flux density sc ale (see e.g., 
Papadoooulos ct al. 2000; Jenness et all 1200 ll) . Given that 
HzRG s appear to be loca ted in submm overdense regions 
(e.g., IStevens et all 120031) flux may have been lost due to 
chopping onto a nearby galaxy. However this very unlikely 
to have affect ed more than a few s ources. 

Following lOmontet al. (2001), Table |2] includes a col- 
umn indicating the quality of the observation. Good quality 
data is indicated by an 'A', whereas poor quality is indi- 
cated by a 'B'. Poor quality reflects bad atmospheric condi- 
tions (e.g. large seeing), short integration time (< 2 sets of 
50 integrations each), or poor internal consistency as shown 
by the K-S test (i.e. the measurements were not consistent, 
but it was impossible to determine which were the outliers. 
In such cases the average of all measurements was used). 

2.2 Potential contamination of the thermal 
submillimetre flux 

Because all our objects are powerful radio galaxies, it is 
important to estimate any synchrotron contribution to the 
observed submm band. We used flux densities from the 
WENSS (325 MHzlRengelink et al.lll997l) . Te xas (365 MHz; 
iDouglas et ail Il99rj) and NVSS (1.4 GHz: ICondon et all 



4 M. Reuland, H. Rottgering, W. van Breugel, and C. De Breuck 



1998) surveys to extrapolate to 350 GHz (850 urn) fre- 
quencies using a power law. For 53W06 9, we extrapolated 
from t he 600 MHz and 1.4 GHz values in IWaddington et alJ 
( 2000). We find that the synchrotron contribution at 850 ii.ni 
is negligible for most galaxies in our sample. Only for some 
objects from A01 (and PKS 1354—17) would this require 
corrections larger than the la uncertainties in the 850 [im 
measurements. 

Synchrotro n spectra often steepen at high frequen- 
cies (e.g., API: lAthreva et al]ll997l lAndreani et aLlEooll 



Sohn et al.ll2003h . and linear extrapolation should be con- 
sidered an upper limit to the synchrotron contribution. A01 
performed parabolic fits to account for the curvature of 
the radio spectrum. Using the midpoint between linear and 
parabolic fits they find corrections larger than 1.5 mjy to 
the 850 urn flux densities in only 6 cases. One could argue 
that parabolic fits are more appropriate, in which case all 
corrections would be negligible. 

Note that the radio measurements reflect the spatially 
integrated flux densities of the sources. The radio cores have 
flatter spectra than the lobes and could dominate at higher 
frequencies. However, they are usually faint and for USS 
sour ces even the radio cor es tend to have steep spectra (e.g., 
API: lAthreva et aljfl997h indicating that contamination by 
the core is likely to be negligible as well. 

Given these uncertainties, extrapolation from the radio 
regime to submm wavelengths is uncertain and is likely to re- 
sult in an overestimate of the non-thermal contribution due 
to steepening of the radio spectrum. This is demonstrated 
for th e case of B3 J233P+3927, for which |Pe Breuck et alJ 
(2PP3) estimate a non-thermal contribution of ~ 1.3 mjy at 
113 GHz but measured a flux density < 0.3 mjy that seems 
to be of thermal origin. Therefore we do not correct for a 
contribution to the submm continuum from the non-thermal 
radio emission, except for a fe w sources discussed b elow. For 
these reasons (and following IWillott et aljliool) we have 
chosen also to use uncorrected fluxes from A01 in the re- 
mainder of this paper. 

There are two exceptions that we exclude from our final 
sample. Following A01, we reject B2 0902+34, as this source 
has a bright flat-spectrum radio core which could dominate 
the submm emission. We also reject the flat spectrum radio 
source PKS 1354—17. It is significantly brighter than any 
of the other sources at 850 [im (5*850 = 20.5 ± 2.6 mjy), 
but linear extrapolation from the radio regime shows that 
the non-thermal contribution at 850 urn could be as large as 
40 mjy and could easily account for all of the submm signal. 

Gravita tional lensing may be important in some cases 
resulting in enhanced submm fluxes. However, 
recen t estimates (e.g., IChaoman et alJl2002bt iDunlop et"aTl 
2002) show that this is limited to a small but significant frac- 
tion (3-5 per cent of sources with 5*850 > lOmJy may have 
been boosted by a factor > 2) and that for most objects 
there is no evidence for strong gravitational lensing. Cor- 
rections for lensing must be made on a case-to-case basis 
and are strongly model dependent. Since these corrections 
are likely to be small, they have not been attempted for the 
present sample. 



2.3 Dust template 

As ha s been noted by many authors (e.g., A01: lHughes et all 
119971) . choosing the dust template is an important step in 
inferring the bolometric far infra-red luminosities (Tfir), 
star formation rates (SFR) and dust masses (Ma). A com- 
plication is that the appropriate dust template may change 
over redshift due to changing dust properties with the evo- 
lutionary states of the galaxies. 

Throughout this paper we adopt single temperature, 
optically thin greybody emission for two sets of emissiv- 
ity index (3 and te mperature T d (see iDunne fc EalesllioOll : 
iDuoac et alJl2003l . for possible concerns) as the functional 
parametrization for thermal dust emission from high red- 
shift sources. We choose j3 = 1.5 and Td = 4 K for com- 
parison with other pap ers (e.g., A01; and see IDunne et alJ 
l200Ct lEales et all 12003ft . and also briefly investigate the ef- 
fects of assuming f3 = 2.0 and Td = 40 K as seems rea- 
sonable f or some hyperlumi nous IR galaxies (Td = 35 K; 
HvLI RGslFarrah et alJl2002ft and z > 4 quasars (T d = 40- 
50K: |Priddev fc McMahor]l200lilWillott et al.lf2002ft . Mea- 
suring the value of /3 for HzRGs directly would require ob- 
servations at many more rest-frame FIR wavelenghts than 
presented here. Increasing /3 or the dust temperature de- 
creases the inferred Lpir of high redshift sources relative to 
lower redshift sources for a given flux density. 

The fraction of absorbed UV /optical light, Ssb, and 
possible departures from the prototype Salpeter ini- 
tial mass function, parametrized with 5imf, are other 
uncertain factors. Generally accepted approximations 



(see e.g. , iPana dopoulos et, alJ l20PPl : lOmont et alJ |2' 
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iDe Breuck et al J 120031 ) for the dust mass, inferred FIR lu- 
minosity and star formation rate are respectively: 



M d = 



5 Q bsTY 



(1 + 2)Kd(^reat)i3(^re S t,T d ) ' 



Tfir = 47rM d 



8nh K d (v) /kT\^ 
~ 2 W \~~h) 



K d {u)B{v,T A )dv = 



r(/3 + 4)C(/3 + 4)M d , 



and 



SFR = <5imf<5sb( Tfir/10 10 L q ) M yr" 1 , 



with Kd(^) oc the frequency dependent mass absorption 
coefficient which modifies the Planck function, B(u,T d ), to 
describe the isothermal greybody emission from dust grains, 
r the Gamma function, ( the Riemann Zeta function, Dl 
the luminosity distance and 5 b s the observed flux density. 
The mass absorption coefficient is poorly constrained (e.g. , 



Chini et alJl98ftlDownes et alll99alPe Breuck et all2003 : 
James et all2002l) and we conform to the intermediate value 
of ft d (375GHz) = 0.15 m 2 kg -1 chosen by A01. 



3 OBSERVATIONAL RESULTS 

24 radio sources were observed. The results of the observa- 
tions, the inferred rest-frame 850 \im luminosities Tsso, total 
far-IR luminosities Lfir, and radio luminosities Lsghz are 
summarized in Tabled 12 of the HzRGs are detected at >3cr 
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significance at 850 \im. The median rms flux density of the 
observations is osso = 1.5 mjy with an interquartile range 
of 0.8 mjy. Only B3 J2330+3927 may have been detected at 
450 \im at a 2a level (S450 = 49.1 ± 17.7 mjy). 

Particularly noteworthy is MG 2141+192. This source 
was detected by A01 at the 4.8cr-level at Ssso = 4.61 ± 
0.96 mjy, using the narrow filterset whereas we observed 
5850 = 2.16 ± 1.10 mjy and S 850 = 2.45 ± 1.58 mjy at two 
separate instances with the wide filter set and did not detect 
the source. Similarly 4C 60.07 has been detected in photom- 
etry mode at Ssso = 11.5 + 1.5, 17.1 + 1.3 and i n a jiggle map 
at 21 .6 ± 1.3 CA01: lPapadopoulos et"aill200Ctlstevens et alJ 
2003). These measurements are consistent to within 2-3<7 
from the mean. However, naively, they could also be inter- 
preted as signs of submm variability. A01 and lWillott et al] 
(2002) also found tentative evidence for variability in the 
submm, but ascribed it to problems with sky subtraction 
for data obtained with the single-element bolometer UKT14 
versus SCUBA. It is hard to see how widespread star forma- 
tion could result in submm variability on the time-scale of 
years. Significant changes in Lfir might be easier to envis- 
age as the result of UV variability often seen in AGN and if 
the submm emission results from quasar heated dust. How- 
ever, even in this scenario changes in the UV are expected 
to average over time in the observed submm regime, unless 
the FIR emitting region is compact. While the typical scale 
sizes for UV emission from the AGN are on parsec scales, the 
minimum extent of the FIR emitting region must be about 
1 kpc to m atch the observed luminosity and dust temper- 
ature (e.e.. ICarilli et alJl200ll) . Submm variability, if real, 
is therefore hard to explain if reprocessing by dust is the 
dominant mechanism for the FIR emission. If, alternatively, 
the FIR emission would be non-thermal emission from the 
AGN, then the emission should be unresolved in contrast to 
the observed extents of a few tens of kpc. Moreover, obser- 
vations indicate that AGN contribute at most 30 per cent 
of the FIR luminosity at wavelengths longer than 50 [i.m, 
at lea st for HyLIRGs iRowan-Robinsonll200Ci iFarrah et alJ 
120021) . 

Possible explanations therefore may be that (i) for dif- 
ferent chopping angles and distances sometimes flux is ac- 
cidentally lost due to companion gal axies in the o ff-beam 
since the fields are overdense iStevens et al]l2003ft . or (ii) 
that sometimes the AGN do in fact contribute close to the 
maximum amount expected (note that the s ubmm flux for 
4C 6 0.07 is centrally concentrated; Fig. 2 in IStevens et alJ 
2003), even though the starbursts still dominate, (iii) Finally 
of course there is still the possibility of pointing errors, un- 
certainties in absolute flux calibration, and differences in at- 
mospheric transparency that could shift the effective band- 
pass by a few GHz, which could make a difference due to 
the very steep slope of the spectrum. 



4 ANALYSIS 

In the following we discuss a sample of 67 radio galaxies 
(46/47 from A01, 23/24 from this paper, two sources were 
observed in both samples). This excludes the flat-spectrum 
sources B2 0902+34 and PKS 1354-17 (see Section|23. For 
the statistical analysis we use the inverse variance weighted 
averages of the measurements of MG 2141+192 (<S850> = 



3.3 ± 0.7mJy) and 4C 60.07 (<S 8 5Q>= 14.4 ± l.OmJy) 
and for MRC 1138—262 we prefer the valu e of Ss50 = 
5.9 + 1.1 mjy obtained bv lStevens et alJ i2003l) over our ob- 
servation under adverse conditions. The median rms flux 
density for this entire sample is osso = l.lmjy with a 
0.24 mjy interquartile range. 

Figure shows that the observed submm flux densities 
and therefore the inferred luminosities (see Table[5J at z > 3 
vary significantly from object to object. For the assumed 
dust template we find a range from Lfir < 4 x 10 12 L© 
for undetected targets to Lfir ~2x 10 13 Lq for detected 
sources. There are several viable scenarios to explain this. 
First, if all the warm dust is heated solely by young stars, 
then Lfir is linked to the SFR, implying that the SFR dif- 
fers significantly between objects. Alternatively, there may 
be a range in produced dust masses as substantial dust pro- 
duction may take more than a billio n years if low-m ass stars 
are the principal contributors fe.g.. lEdm unds 200J). In this 
case the range in Lfir may reflect a range in starburst ages. 
The recent detection of significant am ounts of dust in th e 
local supernova remnant Cassiopeia A iDunne et aljl2003T) . 
indicates that much faster evolving massive stars may be at 
least as important in producing dust. A considerable con- 
tribution from high-redshift supernovae would significantly 
lower the required time-scales for dust production and would 
favor a range in SFRs to explain the range in Lfir. 

Despite the various uncertainties, the result that the far 
infrared luminosities are high is robust. This has important 
implications for the starburst nature of these galaxies: in- 
serting the values for Lfir that we find using either dust 
template confirms that HzRGs are vigorously forming stars 
up to rates of a few 1000 Mq yr _1 and have dust masses of 
a few times 10 8 M© . This is consistent with the notion that 
they are in a critical phase of their formation, forming the 
bulk of their stellar masses. 

4.1 Statistical analysis 

We have performed statistical tests to search for possible 
correlations of the submm with other properties of the ra- 
dio galaxies. Specifically, we investigated whether there are 
correlations with redshift (as reported by A01), radio lumi- 
nosity (as indicator of AGN contribution), largest angular 
size of the radio source (as indicator of age), Lf-band mag- 
nitude (as indicator of stellar mass or star formation rate), 
Lya flux (as a possible indicator of starburst "fuel"), and 
UV polarisation (as indicator of the relative contributions 
of starburst and AGN to the UV continuum). Below, we 
discuss these in more detail. 

Since the sample contains a large fraction (of order 50 
per cent) of non-detections (i.e. upper limits) we have con- 
ducted survival analysis tests similar to A01. Survival anal- 
ysis allows the mixing of detections and upper limits ("cen- 
sored data points" ) in a statistically meaningful way thereby 
preserving as much information as possible. The results of 
the survival tests are summarized in Tables |3| and 0] Table 
^represents the results for the entire sample, whereas Table 
2]only considers the 32 galaxies with redshifts z > 2.5. 

We have defined the subsample with redshift z > 2.5 to 
test for selection effects and remove them from our sample. 
The well known and strong (see A01) radio power-redshift 
relation in flux limited samples is likely to be the origin of 



6 M. Reuland, H. Rottgering, W. van Breugel, and C. De Breuck 



Table 2. Observed 850 \xm and 450 \xm submm flux densities Ssso yxm and S450 yxm with their standard errors for the radio sources in the 
program. The total duration of the observations, N[ n t is given in sets of 50 integrations. 3<r upper limits to the 850 ixm flux are shown for 
sources whose S/N is below 3. Only B3 J2330+3927 may have been detected at 2<r at 450 |jm. Logarithms of inferred rest-frame 850 |om 
luminosities -LssOi far-IR luminosities, LpiR and radio luminosities -L3GHZ are shown for the dust template with (3 = 1.5, T d = 40 K and 
a flat universe with !)m = 0.3, Q\ = 0.7, and Hq = 65 km s" 1 Mpc -1 . 



Source 




z 


iVint 


Ss50 a 


S/N 


Quality 


3<r lim. 


S450 


^850 




AgHz 








x 50 


mJy 






mJy 


mJy 


WHz" 1 sr" 


1 L 


WHz^sr 1 






210 


4+4 


1 Q _1_ 1 Q 

— r.y in 1.0 


1 4 


A 




9 -1- 9Q 




Sli.UO 


94 


MRC1 1 9fS9 


2 


.156 


2 


1 a -|_ q q6 


3.9 


13 




-65 ± 134 


23.26 b 


12. 83 6 


27.15 


WN T1 1 1 ^4-^01 

VV lNJ_Ll_LO-t-OUlU 


2 


oou 


4 


q n _ L 1 q 
o.u H J. .0 








90 A- 1 1 


^9^ ^1 


s 1 9 on 


9^» 


WN 107474-^^^4 


2 


990 


g 


4.8 + 1.1 


4 5 


A 




18 + 15 


23.15 


12.73 


26.22 


VV IN JU^O ItOUUU 


3 


080 


7 


5.9 + 1.6 


3 7 


]3 




-29 ± 22 


23.23 


12.81 


26.27 


.DO J ZOOUTOy^ t 


3 


086 


3 


14.1 i 1.7 C 


8 5 


j\ 




49 + 18 


23.61 


13.19 


26.56 


TN T1 1 1 9—9Q4S 


3 


090 


5 


5.8 + 1.1 


5 1 


j\ 




15 + 9 


23.23 


12.81 


26.66 


MRC0316-257 


3 


.130 


2 


0.6 ± 2.7 


0.2 


B 


<8.8 


5 ± 48 


<23.41 


<12.99 


27.22 


PKS1354-17 


3 


150 


2 


20.5 ± 2.6 d 


8.0 


B 




-47 ± 77 


23.77 




27.71 


WNJ0617+5012 


3 


.153 


6+6 


1.0 ± 0.7 


1.3 


B 


<3.2 


3 ± 16 


<22.96 


<12.55 


26.10 


MRC025 1-273 


3 


.160 


2 


0.6 ± 2.8 


0.2 


A 


<8.9 


-54 ± 91 


<23.41 


<12.99 


27.00 


WNJ1123+3141 


3 


.220 


8 


4.9 ± 1.2 


4.1 


A 




4 ± 14 


23.15 


12.73 


26.60 


WNH1702+6042 


3 


.223 


1 


-0.4 ± 3.6 


-0.1 


B 


<10.8 


-73 ± 91 


<23.49 


<13.07 


26.40 


TNJ0205+2242 


3 


.506 


6 


1.3 ± 1.3 


1.0 


A 


<5.2 


27 ± 23 


<23.17 


<12.75 


26.58 


TNJ0121+1320 


3 


.517 


7+4 


7.5 ± 1.0 


7.6 


A 




4 ± 16 


23.33 


12.91 


26.55 


6CJ1908+722 


3 


.532 


6 


10.8 ± 1.2 C 


9.0 


A 




33 ± 17 


23.49 


13.07 


27.25 


WNJ1911+6342 


3 


.590 


2 


1.3 ± 3.6 


0.4 


B 


<11.9 


-38 ± 50 


<23.53 


<13.11 


26.26 


MG2141+192 


3 


.592 


7+5 


2.3 ± 0.9 b ' e 


2.6 


A 


<5.0 


12 ± 13 


22.96 i > 


<12.55 6 


27.30 


WNJ0346+3039 


3 


.720 


4 


-0.5 ± 1.3 


-0.4 


A 


<3.8 


-5 ± 12 


<23.02 


<12.61 


26.43 


4C60.07 


3 


.791 


5 


11.5 ± 1.5 b - c ' e 


7.6 


A 




10 ± 13 


23.61 6 


13.19 6 


27.15 


TNJ2007-1316 


3 


.830 


5 


5.8 ± 1.5 


4.0 


A 




4 ± 45 


23.21 


12.79 


26.98 


TNJ1338-1942 


4. 


.100 


4+7 


6.9 ± 1.1 


6.2 


A 




-36 ± 32 


23.29 


12.87 


27.05 


TNJ1123-2154 


4. 


.109 


2 


1.5 ± 1.7 


0.9 


A 


<6.7 


-7 ± 11 


<23.27 


<12.85 


26.76 


TNJ0924-2201 


5. 


.190 


8+4 


-0.7 ± 1.1 


-0.7 


A 


<3.2 


-0 ± 26 


<22.94 


<12.53 


27.24 



a This does not include the 10—15 per cent uncertainty in absolute photometric calibration. 

6 For the statistical analysis we use S 8 50 = 5.9±l.lmjy, 5 8 50 = 3.3 + 0.7 mjy and S 8 50 = 14.4+ 1.0 mjy for MRC 1138-262, MG 2141+192 
and 4C 60.07, respectively. Lgso and X/fir were inferred using t hose v alue s. See Sect ion HI for deta ils. 
c Data published in CO imaging studies bv lPap adopoulo s et al ] fcOOrfl and lDe Breuck et alJ i2003ft . 

PKS 1354—17 is likely to be dominated by non-thermal emission (c.f. Section |2~21. 
e Also part of the survey by A01. 



Table 3. Results of survival analyis for all 67 sources. The significance P is the probability of the variables not being correlated according 
to each test. If all three tests yield P < 5 per cent, then the variables are taken to be correlated. If only some yield P < 5 per cent then 
a correlation is regarded as possible, but uncertain. Results arc shown for both a 3<r and 2<r (between brackets) detection limit. 



Variable 




Percentage of 




Significance (P) 




Correlation 


Dependent 


Independent 


Data Censored 


Cox 


Kendall 


Spearman 


Present? 


■Sssoum 


z 


67% (52)% 


0% ( 0%) 


0% ( 0%) 


0% ( 0%) 


YES (YES) 


Ain 


L3 GHz 


0% 


17% 


34% 


32% 


NO 


£3 GHz 


z 


0% 


1% 


2% 


2% 


YES 




£-3 GHz 


67% (52)% 


1% ( 1%) 


7% ( 8%) 


10% ( 8%) 


MAYBE (MAYBE) 


L&50\xm 


Z 


67% (52)% 


0% ( 0%) 


0% ( 0%) 


0% ( 0%) 


YES (YES) 




a 


67% (52)% 


1% ( 1%) 


1% ( 0%) 


1% ( 0%) 


YES (YES) 


£850um 


Ain 


67% (52)% 


56% (22%) 


79% (89%) 


96% (78%) 


NO (NO) 


a 


L3 GHz 


0% 


36% 


20% 


19% 


NO 


a 


Z 


0% 


0% 


0% 


0% 


YES 



many apparent correlations with redshift (see Table [3J. An- 
other example of selection effects is the correlation found be- 
tween Lg5o and radio spectral index a (Table |3J. It reflects 
the tight correlations between spectral index and redshift 
(part of our search criteria for HzRGs) and between Lsso 
and redshift. In the z > 2.5 sample this correlation almost 
disappears even though there is a large range in both Lsso 
and a, confirming that it was spurious. An additional advan- 



tage of this subsample is the homogeneity of the supporting 
data (Lya fluxes and A'-band magnitudes are understand- 
ably sparse for lower redshift sources). 

Below, we describe the details of the survival analysis 
and the criteria used to determine whether a correlation is 
present. In Section |5J we discuss the correlations individu- 
ally. 
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Table 4. Similar to Tabled Results of survival analysis for 32 sources at redshifts z > 2.5. 



Variable 




Percentage of 




Significance (P) 




Correlation 


Dependent 


Independent 


Data Censored 


Cox a 


Kendall 


Spearman 


Present? 




z 


53% (34%) 


95% (81%) 


17% (25%) 


34% (38%) 


NO (NO) 


Ain 


£3 GHz 


0% 


66% 


8% 


9% 


NO 


^3 GHz 


Z 


0% 


7% 


4% 


4% 


MAYBE 




£3 GHz 


53% (34%) 


21% (33%) 


33% (48%) 


46% (58%) 


NO (NO) 




Z 


53% (34%) 


88% (75%) 


20% (34%) 


35% (52%) 


NO (NO) 




a 


53% (34%) 


88% (99%) 


81% (69%) 


99% (71%) 


NO (NO) 


Ls50\Ma 


Ain 


53% (34%) 


13% ( 6%) 


16% ( 7%) 


13% ( 6%) 


NO (NO) 


Ls50\Ma 




59% (41%) 




86% (91%) 


4% ( 4%) 


MAYBE (MAYBE) 


i850|jm 


Kmag 


59% (44%) 




47% (33%) 


88% (72%) 


NO (NO) 


a 


£3 GHz 


0% 


56% 


21% 


11% 


NO 


a 


Z 


0% 


2% 


10% 


8% 


MAYBE 



a Cox's proportional hazard model only allows censoring in the dependent variable, therefore this test could not be applied 
with the Lyo and if -band data. 



Survival analysis 
The application of survival analysis methods to 
astronomical data has been described i n deta il by 
iFeieelson fc Nelson! l|l985h and llsobe et alJ Jl98ftt . We 
have m ade use of t he routines in the STSDAS package of 
IRAF 2 jTodvlll993l), which we re modelled after the ASURV 
package ijLavallev et alJ 1992). Cox's proportional hazard 
model, the generalized Spearman's rank order correlation 
coefficient, and the generalized Kendall's tau correlation 
coefficient test the null hypothesis that no correlation 
is present in the sample. We adopt the convention that 
two variables are correlated if the chance P of the null 
hypothesis being true is smaller than 5 per cent. While these 
tests should give similar results they have different specific 
limitations and strong points. Therefore, a correlation is 
considered to be reliable if all tests yield P < 5 per cent, 
and possible but unconfirmed if only some indicate P < 5 
per cent. Note that the Cox test only allows censoring in 
the dependent variable and that the generalized Spearman's 
Rho routine is not reliable for small data sets (N < 30). 
SCUBA occasionally yields negative flux densities. We 
defined the na upperlimit for those cases to be nx the 
rms flux density. Similarly, the no upperlimit for a positive 
signal is defined as S + nx the rms flux density, with S the 
observed signal. 




Redshift 



Figure 1. Observed 850 um flux density and 3cr upper limits 
versus redshift of all 67 radio galaxies discussed in this paper. 
The size of the arrows and errorbars correspond to lc rms. 



5.1 Redshift dependent submillimetre properties 



5 CORRELATIONS BETWEEN PARAMETERS 

Tables |21 and 2] represent the results of the survival analysis 
as described above. We will now discuss the motivation for 
each test and the implications of the correlations (or lack 
thereof) found in order of importance. 



2 IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



5.1.1 Flux density and relative detection fraction 

Table |3 shows that the observed submm flux densities Ssso 
are strongly correlated with redshift. This is reflected in Fig- 
ure and [5] Figure shows the observed 850 \im flux den- 
sities against redshift. Figure |2] shows the number of radio 
galaxies that have been observed in a particular redshift bin 
and the number of those which have been detected given a 
2(7 or 3<r detection criterion. The success rate of our z > 3 
program is truly remarkable, and corroborates the conclu- 
sion of A01 that the detection fraction of ~ 50-67 per cent 
at z > 2.5 is significantly different from the detection frac- 
tion of ~ 15 per cent at z < 2.5. Removing the 4 brightest 
sources from the sample does not destroy this relation. 
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Figure 2. Histograms of 850 um SCUBA detections with S/N > 
2 (left) and S/N > 3 (right) versus the total number of observed 
radi o galaxies (grey) as a function of redshift. The detections 
from I Archibald et alJ l200ll> are shown as dashed, the detections 
from our program are shown in black. At z > 2.5 ~ 50-67 per 
cent of the galaxies are detected, as opposed to ~ 15 per cent at 
z < 2.5. The bins have a width of unit redshift and are centered 
at z = 1,2,3,4,5. 



5.1.2 Investigation of the difference between detections 
and non- detections 

Given a 3a detection criterion the average submm flux den- 
sity of the 22 detected radio sources and 45 non-detections 
are <Sgso,^3<?> = 6.60 ± 0.71 mjy and <S850,<3 CT >= 0.81 ± 
0.20 mjy respectively, while the average for the entire sam- 
ple is <S , 850,sam P ic>= 2.71 ± 0.43 mjy. For a 2a detec- 
tion criterion the average submm flux density of the 30 
detected galaxies and 37 non-detections are <Ss5o,^2<y>~ 
5.29 ± 0.60 mjy and <S 8 50,<2<r>= 0.34 ± 0.19mJy respec- 
tively. 

Figure [3] shows cumulative redshift distributions £(z) 
for various subsets of the radio galaxy sample together with 
spectroscopic redshifts for other submm samples from the 
literature. As was suspected from Figure [5] the radio galax- 
ies detected in the submm follow a distribution that differs 
significantly from both the undetected sources and the par- 
ent sample. The K-S test shows with P > 99 per cent con- 
fidence that using a 2a or 3a detection criterion picks out 
the same population (both for the detections and for the 
non-detections), whereas the chance that the detections and 
non-detections are distributed similarly is P < 1 per cent. 

Does the high median redshift (z = 3.1) of the detected 
sources (the parent sample has z = 2.5) purely reflect the 
strong negative if-correction or does it reflect a change in 
-i/FiR? Figure shows that the median redshift of the de- 
tections is highe r than the median red shift (z = 2.4) for 
SCUBA sources JChapman et al]l2003l) . while the redshift 
distribution of the parent sample follows the SCUBA popu- 
lation closely. This argues in favor of a different I/fir for de- 
tections and non-detections. What then causes the difference 
in redshift distribution between detected HzRGs and the 
submm population? There are at least two possible explana- 
tions. It may indicate th at HzRGs are more ma ssive, at the 
centers of protoclusters dVenemans et al]l2003l) . and there- 
fore undergo a faster evolution than 'normal' submm sources 
and finish the bulk of their formation process early. Alter- 
natively, it may reflect a selection effect as the requirement 
of a faint radio counterpart prior to spectroscopic follow-up 
for the submm sample selects against high redshift galaxies 
JChapman et alJl2003t . 




2 4 
Redshift 



Figure 3. Cumulative histograms T,N(z) for redshifts of submm 
sources from the literature and the 67 radio galaxies (solid line) 
discussed in this paper. The short-dashed and the dotted line in- 
dicate submm detections of radio galaxies versus non-detections 
given a 3cr detection criterion, respectively. The detections have 
a higher median redshift (z = 3.1) than the parent sample (z = 
2.5). The two long-dashed lines represent £7V(z) of the 10 submm 
source s with spectroscopic redshifts found by |£!haDrmme^il] 
( 2003) and using 9 add itional spectroscopic redshifts ( Ivisoiiet^L| 



19981 l200rilDev et alll99fltlSoucail et alJltWfltlEales et al.l2000 
all2002UChapman et al.l2002blcUAretxaga et al.l2003 
il*Jl2003l) . 



Lcdlow et 
iFraver et al 



5.1.3 The increase of submm luminosity with redshift 

Figure 0] shows the inferred 850 [im luminosity Lsbo as a 
function of redshift. For each redshift bin Lsso has been es- 
timated using the Kaplan-Maier estimator, a survival analy- 
sis technique which tries to estimate the true distribution of 
the underlying population by incorporating both upperlim- 
its and detections. We have computed the luminosities both 
using a dust template with (5 = 1.5 and (3 = 2.0. Adding our 
data to the sample of A01 confirms that HzRGs have higher 
submm luminosities than lower redshift sources. However, 
there are significant differences with the findings from A01: 
(i) 1/850 at low redshifts is higher than inferred by A01. This 
is because we have chosen, not to correct for possible syn- 
chrotron contamination (see Section 2.3) (ii) There is weak 
evidence for a turnover or leveling off at z > 4 (based on the 
lower flux densities of the five galaxies with z > 4 compared 
to redshifts 3 < z < 4; see Fig. Ql. This explains why there 
is no strong correlation between Lsso and z in the high red- 
shift subsample (see Table [IJ. (iii) A01 remarked that the 
increase in Lg.50 becomes less pronounced if one assumes a 
dust template with (3 = 2. Because Lsso at low redshifts is 
higher than in A01, this effect is even stronger in our sample 
and we find that Z/ 8 50 may be rather constant. The sources 
would still be highly submm luminous, implying huge dust 
masses and star formation rates. 
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Figure 4. The average rest-frame 850 |om luminosity versus red- 
shift for 67 radio galaxies as estimated using the Kaplan-Maier 
estimator for each redshift bin. The data points indicated with 
open squares and circles were calculated assuming the same dust 
template as A01 (/3 = 1.5, T d = 40 K) for 3<r and 2cr detection 
criteria respectively. The solid symbols are similar but assume 
j3 = 2.0 instead. The data points at z = 5 are 3<r upperlimits 
based on the non-detection of TN J0924-2201. 



5.2 The connection between submm and radio 
luminosity 

Radio galaxies host luminous AGNs and for favorable ge- 
ometries (e.g., dusty warped disks) heating of dust by 
only a small fraction 20 per cent) of their UV radia- 
tion could easily ex plain typical far-IR luminosities (e.g., 
ISanders et all [1989). Therefore, an important question is: 
do the AGN heat the large-scale dust significantly, i.e. are 
iradio and Lfir correlated? 

To answer this question we have estimated the rest- 
frame radio power at 3 GHz, Z/3GHz- We have chosen this 
frequency because for the redshift range 1.1 < z < 7.2 this 
requires only interpolation between between 365 MHz and 
1.4 Ghz. Of course this assumes that the UV e mission of the 
central source relates to its ra dio output fe.g.. lWillott et alJ 
ll999l:lDe Breuck et al.ll2000b|) . 

While Table [3] shows that there is a possible correlation 
between Lsso and Z/3GHz over the entire redshift range of 
our sample, this is probably an effect of the strong redshift 
dependence of both Lg$o and Z/3ghz- The tests represented 
in Table 2] only take radio galaxies with z > 2.5 into ac- 
count and even though this subsample contains almost all 
submm detections the correlation disappears. The scatter- 
plot in Figure |S] is a graphical representation of the same 
result. We don't find evidence for dust heat ing by AGN, 
in accordance with earlier findings (e.g., A0 1;|Will ott et alJ 
120021: lAndreani et aljl20p3> but see lHaas et alJ fcOOSl) for a 
contrasting view based on ISO observations of local (most 
have 0.1 < z < 1.0) hyperluminous quasars. This analysis is 
likely to be an oversimplification because the typical time- 
scales for starburst and radio activity differ by an order of 
magnitude. 



22.5 L- 1 
25.5 
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Figure 5. Submillimetre luminosity against radio power L3GHZ 
for 32 radio galaxies with z > 2.5. No strong correlation is appar- 
ent, indicating that heating of the dust by AGN is not a dominant 
process for radio galaxies. Symbols as in Figure 111 



5.3 An anti-correlation between submm flux and 
UV polarisation 

Optical polarimetry studies can be used to search for di- 
rect signs of starburst activity and determine the relative 
contr ibutions of AGN and starburst to t he U V continuum 
(e.g., iDev et al.l|l997l : IVernet et ail 1200 lfl . If the starburst 
dominates the UV/optical light, then the scattered (i.e. po- 
larized) AGN light is diluted and one expects a low polarisa- 
tion fraction of the UV continuum and perhaps a correlation 
with the observed submm flux densities. 

Pol arimetry results exist for 13 of the galaxies in our 
sample (IDev et al.lll997t IVernet et ail 1200 ll De Breuck et 
al. in preparation) . A summary of the polarisation fractions 
and corresponding 850 u.m flux densities is given in Table [5] 
We have plotted the observed submm flux density against 
UV continuum polarisation fraction in Figure |S] No sources 
with high UV polarisation are detected in the submm. This 
is supported formally by survival analysis: Kendall Tau gives 
a probability P — 3 per cent that there is no correlation. 

This result is consistent with the view that for HzRGs 
star formation and submm emission are closely linked 
while any AGN contribution to the FIR is negligible. 
iTadhunter et alJ (I2002T) also found some evidence for this. 
Two out of the three objects which are starburst dominated 
in their sample of 2-Jy radio galaxies are also extremely 
FIR luminous. Further support for low polarisation in star- 
bursting systems co mes from the luminous s ubmm source 
SMM J02399-0136 iVernet fc CimattilioOlh which shows 
only moderate (P ~ 5 per cent) polarisation. 

The weakness of the submm luminosity in highly polar- 
ized sources together with the tentative correlation between 
Lycv and submm flux (see next section) is consistent with 
the anti-correlation between Lya l umino sity and UV polar- 
isation reported bv IVernet et al] (l200lT) . They argue that 
Lya photons are resonantly destructed by dust and that 
scattering by the same dust results in a higher fractional 
polarisation. This explanation, however, seems at odds with 
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Table 5. UV/optical polarisation fractions and submm fluxes 
for all radio galaxies for which both ha ve been observed . 
UV/optical polarisa tion data are taken from lDev et alJ ll997l) , 
IVernet et alj 1200 ll) . and De Breuck et al. in preparation, the 
submm data are from A01 and this paper. 



Source z P(%) S 850 (mjy) 
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< 1 


4.1 
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2.483 


15.3 
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3.130 


< 4 
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10.8 
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Figure 6. 850 um flux 5850 against observed UV continuum po- 
larisation fraction. None of the highly polarized sources are de- 
tected in the submm, indicating that low polarisation (indicative 
of a starburst) and high LpiR both trace star formation. 

the lower submm flux densities in highly polarized sources 
unless there would be significant amounts of dust that is too 
hot to be detected in the submm. 



5.4 Submillimetre and Lya flux 

Since many HzRGs are embedded in giant Lya haloes (e.g. , 
iMcCarthvllTiffli Ivan Oiik et all 11999 : iReuland et al.ll2003fi 

it would be interesting to see if one can relate the amount of 
available star formation 'fuel', as estimated roughly by the 
spatial size, or luminosity, of these emission line haloes and 
the amount of HI as probed by Lya absorption to the ob- 
served rest-frame FIR luminosities. Alternatively, one might 
expect a strong anti-correlation due to destruction of Lya 
emission by dust. 

Many HzRGs are identified based on their Lya line 



and good data on the observed Lya fluxes is available in 
most cases. There is a strong selection effect, however, since 
sources with low Lya flux densities are less likely to be recog- 
nized as HzRGs. Interestingly, Tableland Figure[7|indicate 
that a correlation may be present. While this relation is ten- 
tative only, it offers perspectives for future programs. Better 
data should make it possible to investigate whether there is 
a possible relationship between the dust content and Lya 
to C IV or N v emission line ratio, as indicators of metalicity 
(c.f. IVernet et al]l200lfl . T his would be especial l y inte rest- 
ing in the l ight o f findings bv lHamann fc Ferlandl <ll999l) and 
iFan et ail (j200lT l that QSOs show (super)solar abundances 
out to the highest redshifts and must already have under- 
gone significant star formation. A comparison with the more 
readily studied HzRG host galaxies would be crucial to bet- 
ter understand this result. 

The detection of both Lya and dust in many HzRGs 
is intriguing because already a small amount of dust can 
extinguish Lya radiation efficiently. The detect ion of Lya 
in the submm sources of IChapman et ail J2003T) is equally 
surprising. One explanation may be that Lya emission and 
dust are located predominantly in different locations (e.g., a 
merger between one dusty and one less dusty component). 
Some evidence for this comes from 4C 60.07 at z = 3.8 
in which the spatially re solved dust continuum emission 
iPapadopoulos et alj|2 000|) is anti-corr e lated with the Lyq 
emission (Fig. 5 in lReuland et aDl2003ft . ISmail et al.l (120031) 
also report evidence for an offset between far-IR and UV 
emitting regions in the z = 2.38 starburst galaxy N2 850.4. 

5.5 isso and linear size 

IWillott etall (2002) found an anti-correlation between lin- 
ear size and 850 urn luminosity for their sample of z ~ 1.5 
radio loud quasars. This effect was attributed to a possible 
relation between the jet-triggering event and a short-lived 
starburst or quasar-heated dust in QSOs. No such corre- 
lation was found for a matched subset of the A01 sample 
centered at the same redshift. We have briefly investigated 
whether such a correlation exists in the present sample. 
Figure |H| shows the observed submm flux plotted against 
the projected linear size. No correction for possible differ- 
ences in inclination angle was made because such correc- 
tions are expected to be relatively small among radio galax- 
ies (even when comparing radio galaxies with quasars pro- 
jection effects result in only a factor ~ 1.6 difference in 
proje cted linear size for (9trans = 53 ° ; IWillott et alJ l2000l . 
2002). Our analysis does not support such a strong rela- 
tion and if one interprets linear size as an estimate of the 
age of the radio source (reasonable if one assumes that the 
environ ments and jet powers do not change significantly; 
sec e.g., iBlundell et~aT] ll99Sl^ . then this would be consis- 
tent with the scenario sketched by Willott et al. that the 
correlation found for QSOs is indicative of short starbursts 
while the starbursts in HzRGs have longer time-scales and 
could be forming the bulk of the stellar populations. De- 
spite the absence of a clear correlation, some relation be- 
tween age of radio source and starburst might have been 
expected if the blac k-hole and stellar bul ge grow in a sym- 
biotic fashion (e.g., IWilliams et al.lll999l) . However, in re- 
alistic scenarios many complicating factors can be envis- 
aged such as the significantly different time-scales involved. 
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Figure 7. Sgso um against observed Lya line flux (left) and inferred luminosities Lsso against L^ ya (right). Tentative evidence for a 
correlation is apparent. The galaxies brightest in Lya seem to have higher submm flux densities, possibly indicating that galaxies with 
richer gaseous environments (as traced by Lya) also host more massive starbursts. Symbols as in Figure HI 




D„„ [kpc] 

Figure 8. 850 \xm luminosity Lsso against the radio source pro- 
jected linear size D\j n . There is no evidence for a strong anti- 
correlation as was found for radio loud quasars by Will ott et alj 
(2002). Symbols as in Figure [l] 

Growth of the radio source may fi rst induce star formation 
by compressing mol ecular clouds feegelman fc CiofBlll989t 
iBicknell et alJl200ff) and may then actually signal the end 
of the starbur st by clearing o ut all the cool gas from the 
central region llRawlingsll200l) . 

5.6 Submillimetre and near-IR emission 

As discussed by llsaak et alJ (|2002h . for QSOs one might 
naively expect a correlation between the submm and op- 
tical flux densities. This correlation might arise regardless 
of whether stars or a buried AGN are responsible for heat- 
ing of the dust since black h ole mass scales with stellar bulge 
mass. iPriddev et alJ ll200otl note that the picture is likely to 



be more complex and that any correlation might be smeared 
out due to differences in relative timing between AGN fuel- 
ing and starburst, or varying dust-torus geometries. Radio 
galaxies might be better suited for such studies, sinc e accord- 
ing t o orientation based unification schemes (e.g., iBarthell 
1989) the AGN are obscured by a natural coronograph and 
rest-frame 13-band luminosity therefore is a fairly reliable 
measure of the mass of the stellar population or starburst 
activity. For this reason we have correlated the Zf-band 
magnitudes (samples rest-frame B,V for z > 3) with the 
submm flux densities. As can be seen from Tableland Fig- 
ure HI] no such correlation is apparent. This remains true, if 
we apply a A'-correction to the near-IR magnitudes. How- 
ever, the dataset is rather inhomogeneous and more sensitive 
data in both the submm and near-IR are required to truly 
rule out any correlation. Sensitive multi-band observations 
could be used to (i) construct a redshift independent esti- 
mate of the line-free rest-frame optical continuum for com- 
parison with I/fir, ([Oni] and Ha can sometimes dominat e 
the Zf-band e.g., 4C 30.36, 4C 39.37: lEgami et al.l 12003^ 1. 



and (ii) check whether FIR emission is stron; 
with higher intrinsic reddening (following e.g. , 


?er for sources 


ICalzettil997t 


Adclbcreer & Stcidel 2000; Seibert et al. 


2002). 



5.7 Summary 

From our analysis we consider the following results the most 
interesting: 

(i) Ssso and z are strongly correlated. Whether this is true 
also for Z/fir and z depends on the assumed dust template. 

(ii) -1/850 and radio power L3GHZ do not correlate strongly, 
indicating that AGN do not dominate the submm emission 
(either through synchrotron contamination or dust heating 
by the UV continuum). 

(iii) 5850 anu Lya appear weakly correlated. While this 
is not very convincing, the interesting result is the absence 
of a strong anti-correlation as expected naively from the 
destruction of Lya emission by dust. 
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Figure 10. The 850 um detection (left S/N > 2; right S/N > 3) fractions of radio galaxies as function of redshift (solid line). The 
fractions rise up to a redshift of z = 4. For comparison we also plot the detection fraction for observation s of radio quiet QSOs at z = 1 — 
5, and radio loud quasars at z ~ 1.5 (after matching the sensitivity to 3mjy rms) at 850 um (dashed line:[p"rid deY et al.|2003||]^aak et alj 
l2002UWillott et alJl2002^ and 1.25 mm (dotted line: lQmont et al J20011 120031 ICarilli et all200lllpetric et all2003t iBertoldi et al.l2003h . 

The detection rates for the QSOs appear to follow the same trend as for the radio galaxies, but the QSO datapoint at z = 5 seems 
indicative of a turnover (the 2 = 5 point for the HzRGs is based on only 1 galaxy, but taking all z > 4 HzRGs yields a similar result). 
The datapoint at z = 6 could indicate a detection fraction larger even than what is observed for z = 4. This might reflect, the extreme 
properties of those sources. 
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Figure 9. Ssso um against observed i\-band magnitude. No cor- 
relation is apparent, and no relation between Lpm and stellar 
mass/optical star formation rate of the host galaxy can be estab- 
lished. Symbols as in Figure HI 



(iv) Ssso and UV polarisation fraction appear anti- 
correlated. This is expected if the UV continuum of the 
buried AGN does not contribute significantly to the dust 
heating. 



6 A COMPARISON OF RADIO GALAXIES 
WITH QSOS 

The optically thin (sub-)mm emission of starbursts should 
be largely independent of viewing an gle. Orientatio n based 
unification schemes for AGN (e.g., iBarthell Il989l) there- 
fore suggest that we can combine observations of radio 
galaxies with those of QSOs, allowing us to study the 
evolution of star formation in their host galaxies over a 
larger redshift range. Extensive surveys of high-redshift 



(mostly radio 


-quiet) QSOs have been published recently 


(Carilli et alJ 12001: Omont et alJ 1200 ll 120031: llsaak et al. 


2002 
120031 


iPriddev et alJl2003l iPetric et alJl2003l iBertoldi et al. 


. We have combined their results and compare those 



with our sample. 

Figure [THl shows the detection fraction as a function of 
redshift for the QSO and radio galaxy samples. In many 
redshift bins only a small number of detections and non- 
detections is responsible for estimating the success rates for 
the present effective flux limits. We have estimated the range 
over which the 'true success rate' could vary such that there 
would be a 68 per cent chance of measuring the observed 
detection fraction. This range is indicated by the errorbars. 
Figure ITTI shows the average (sub-)mm flux density for the 
various samples as a function of redshift. 

Reassuringly, the trends for the QSOs and radio galax- 
ies shown in these figures mimic each other. The average 
observed (sub-)mm flux density rises to redshifts of z ~ 4. 
For higher redshifts there is some evidence for a decline (see 
the discussion in Section [5.1.31 and the z = 5 data point of 
the MAMBO QSO obervations) . However, the observations 
of z > 4 QSOs would also be consistent with a fairly con- 
stant or even a rise in the detection fraction and average 
FIR luminosity out to the highest redshifts. This depends 
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Figure 11. The average measured (sub-)mm flux density binned 
in redshift has been plotted for QSOs and radio galaxies. The 
squares represent observations of radio galaxies at 850 [am, the 
triangles QSOs at 850 um, while the circles represent observa- 
tions of QSOs at 1250 um multiplied by a factor of 2.5 for easier 
comparison. The datapoint atzfs 1.5 represents the radio loud 
quasars from lwillott et alj 120021) . The trends with redshift are 
similar to what is seen in Figure Hoi and confirm the steady rise 
from z = 1 to z = 4 for the galaxies. 

rather critically on the z ~ 6 QSOs, which are likely to be 
atypical even for QSOs. 

The average flux density and detection rate of the QSOs 
is comparable to that of the HzRGs, even though the HzRG 
surveys are more sensitive than the QSO surveys by a factor 
of approximately 3. This, coupled with the assumption that 
quasars are on average more luminous than radio galaxies 
in A GN unification schemes based on receding torus models 
fe.g.. ISimpsonl l2003'l . may be viewed as further support for 
an AGN related compo nent in QSOs, similar to what was 
found for radio quasars IIAndreani et al. 2002; Willoti^^ii 
2002). Note that the quasars discussed by IWillott, et al. 
(2002) are exceptionally bright compared to the mostly ra- 
dio quiet QSOs represented in Figs. 1101111 

Due to a lack of carefully matched surveys fundamental 
differences between type I and type fl AGN and radio-loud 
versus radio-quiet targets cannot be properly investigated 
yet, although slowly progress is being made. 



7 SUMMARY 

We have presented SCUBA observations of 24 radio galaxies 
and compared those with earlier results of 47 radio galax- 
ies from the survey by A01. We confirm that HzRGs are 
massive forming galaxies, forming stars up to rates of a few 
thousand Mq yr _1 and that there is no strong evidence for 
a correlation with radio power. Further evidence for a pre- 
dominantly starburst nature of the far-IR emission comes 
from the striking anti-correlation between submm flux den- 
sity and UV polarisation (Fig. 0. 

In agreement with A01 we find that submm detection 
rate appears to be primarily a function of redshift. If this is 



interpreted as being due to a change in the intrinsic far-IR 
luminosity, it would be consistent with a scenario in which 
the bulk of the stellar population of radio galaxies forms 
rapidly around redshifts of z = 3 — 5 after w hich they are 
more passively evolving (c.f. iBest et al"lll99ct) . We also find 
that the median redshift of the HzRGs with SCUBA detec- 
tions (z — 3.1) is higher th an the median redshif t of the 
submm population (z = 2.4: IChapman et alJl2003Tl . In the 
current picture of hierarchical galaxy formation, this could 
be interpreted as that HzRGs are more massive galaxies, 
which are then thought to begin their collapse at earlier cos- 
mic times and evolve faster and finish the bulk of their for- 
mation process earlier. Alternatively, it could indicate that 
higher redshift submm sources are being missed due to the 
requirement of a radio counterpart prior to spectroscopic 
follow-up. 

HzRGs have accurately determined redshifts and host 
identifications and are thought to be the m ost massive 
galaxies at any epoch iDe Breuck et alJ |£o02) . Therefore, 
HzRGs are a key population for studies of galaxy forma- 
tion in the early universe, allowing detailed follow-up mm- 
interferometry observations to study their dust and gas con- 
tent. Currently, they offer the best way to obtain reliable dy- 
namic masses for a significant number of massive high red- 
shift galaxies. These HzRGs would be especially suited to 
constrain any evolution in galaxy mass with redshift, study 
changes in evolutionary status, gas mass, and the starburst- 
AGN connection. If all of these turn out to have masses 
larger than 10 11 Mg, then this could have important con- 
sequences (depending on rather uncertain correction factors 
for the fraction of similar galaxies for which the black hole 
is dormant) for our understanding of galaxy formation, be- 
cause only few such massive gala xies are expected at such 
high redshifts (iGenzel et alJl2003l) . 
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